Counting the negative eigenvalues of the thermalon in three dimensions 
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Some years ago it was shown that the cosmological constant may be reduced by thermal produc- 
tion of membranes that, after nucleation, collapse into a black hole. The probability of the process 
' was calculated in the leading semiclassical approximation by studying an associated Euclidean con- 

figuration called the thermalon. Here we investigate the thermalon in three spacetime dimensions, 
describing the nucleation of closed strings that collapse into point particle singularities. In this 
context we may analyze the one-loop structure without the well known problems brought in by the 
propagating gravitational degrees of freedom. We found that the coupling to gravity may increase 
Q ■ the number of negative eigenvalues of the operator. 
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I. INTRODUCTION 

Currently some of the major challenges of theoretical physics are related to the cosmological constant. First, the 
O |' infamous cosmological constant ■problem j _va which the observed cosmological constant is smaller in at least 60 orders 
of magnitude than the estimated value [l| . Second, there is the question of why the cosmological constant was big at 
i ■ the early stages of cosmic history, producing inflation, and how it relaxed to a very tiny number transferring its energy 
into normal matter. Third, the universe seems to be accelerating [2|, which is normally explained by the presence of 
a positive cosmological constant or, more generally, some kind of (unknown) dark energy whose nature has escaped 
k>" \ our present knowledge. 

In this work we deal with the first two problems. Specifically, we continue the study started in [3] on the decay of 
the cosmological constant by thermal nucleation of bubbles of true vacuum. This process may be described by the 
, so called thermalon, an Euclidean, time independent, solution of a charged membrane coupled to gravity and, in D 
dimensions, a (D — l)-form electromagnetic potential. The nucleation rate of this process is of the form 



T = Ae-», (1) 

o . 

where B — Iback — Ithermt Itherm is the action of the system evaluated on the thermalon, and Iback is the action 
evaluated in the false vacuum background. Here we use the same convention used in in which the Euclidean 
path integral includes e 1 . The prefactor A corresponds to the one loop contribution to the decay rate, namely, it is 
determined by the determinant of the operator O that arises from the second order perturbations of the action 

S (2) Itherm = ~\ J d^xrbO^. (2) 

Here ip stands for the independent degrees of freedom of the action. 

Unfortunately, to compute the above expression is in general extremely difficult. First, one must face the inherent 
difficulties in computing functional determinants. Second, gravity is not renormalizablc for d > 3. This second 
problem may be avoided by analyzing the thermalon solution in three spacetime dimensions, which is what we do in 
what follows. We expect that this toy model may shed some light on the problem of the negative eigenvalues around 
the thermalon. Depending on the classical solution considered the operator O may have any number of negative 
eigenvalues. On the other hand, on some general grounds, Coleman proved [4j that a classical solution which would 
yield an operator O with more than a single negative value "has nothing to do with decay rates" . Nonetheless, that 
demonstration manifestly excludes de Sitter spaces, and therefore there is no objection in principle to extend the 
analysis to operators with a larger number of negative eigenvalues. There is some evidence, however, that Coleman's 
result may be extended to a more general family of theories, including, in particular, field theory in dS backgrounds 
(see [Bj]). We thus expect that if the thermalon is going to describe a decay process, it should have one and only one 
negative eigenvalue. Unfortunately, it is known that in some cases with fixed background the number of negative 
eigenvalues should be at least five[6j. The immediate question is whether the introduction of the gravitational degrees 



2 



of freedom may reduce the number of negative eigenvalues. We explore this possibility in this work and conclude that, 
under some general assumptions, gravity actually increases the number of negative eigenvalues. 

The plan of the article is as follows. In section [ill we review the thermalon general scheme and describe the three- 
dimensional case. Next, in section HTT1 we analyze the eigenvalues of the operator describing second order perturbations 
neglecting the back reaction of the geometry. Then, in section IIVI we analyze the decay rate to first loop including 
the gravitational degrees of freedom. Finally we draw out some conclusions. 

II. THE THERMALON IN THREE DIMENSIONS 

We are interested in the following action, 

Ie = T- I d3x Vd(R - 2A) + 4 / d 3 x^F^ p F^» - f *F AA-q f A- ,i f Vhd 2 x, (3) 

where M. is the bulk and £ the worldsheet swept by a closed string of charge q and tension fi. The electromagnetic 
potential A is a 2-form and F = dA. The third integral is along the boundary of the manifold, and it is necessary to 
have as boundary condition the fixing of F. Finally, we have defined k = 8ttG and units in which AnG — 1 such that 
k = 2. 

The thermalon discussed in this work is a static, spherically symmetric extreme of action above (j3J), therefore we 
will consider metrics of the form 

ds 2 ± = h 2 ± {r)di 2 + f± 2 (r)dr 2 + r 2 d(/) 2 , (4) 

where < t < (3, < <j> < 2ir. The subindices ± indicate the two regions of spacetime separated by the worldsheet 
of the string. The worldsheet may be parameterized such that r = R(t,4>). In this case, however, we are searching 
for a static and axisymmetric solution, thus R(t,<p) = Rq (constant). Since F is a three form in three dimensions 
necessarily F^ p = ^/gEe^ upi where £ is a scalar, allowing to operate in term of the scalar E. For instance the 
electromagnetic equations yield 

E+ - £_ = q. (5) 

On the other hand, the Einstein equations have a matter term coming from the electric field which just adds a piece 
to the cosmological constant. At each region we have 
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A± = X + AirGEi = w . (6) 
*± 

To simplify the calculations, we are going to assume A = 0, such that the cosmological constant come entirely from 
the electromagnetic field. 

Considering the restrictions above the Einstein equations imply that at each side of the string the geometry is 
locally de Sitter. The (+)-region will be assumed to be free de Sitter space. The (-)-region will be the geometry of a 
particle singularity [7] placed at r = 0, 

ds 2 ± = C 2 ± f 2 ± {r)dt 2 + f ± 2 (r)dr 2 + r 2 d^ 2 , (7) 

with 

#(r) = l-£.£(r) = l-^-£. (8) 

Here M is the mass of the particle and C± are arbitrary constants, whose value is determined by the Israel matching 
conditions [1], which in this case require goo to be continuous across the string, hence C_ = f + {R)/f-{R)C + . 

Finally, one can always set the time coordinate, t, such that C+ = 1 and therefore the killing vector <9 t be unitary 
at r = 0. 

The remaining Israel matching conditions give 

f+(R)-f-(R)=2nR, (9) 
while the fact that we are searching for a static solution (R. = R = 0) yields 

f'(R)-f'_(R)=2v. (10) 
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The two unknown quantities above, Rq and M, can be determined by solving equations ([9]) and (fTO]) . The radius 
of the thermalon is given by 



where x = 7 + yl + 7 2 , 7 = Z+ ((2/i) 2 — a 2 ) /4/x and a 2 = — \/l 2 _. Analogously, the mass, M, is given by 

M=*0!*±. (12) 

x 

The geometry in the (— )-region is that of the particle singularity introduced in 0,0- To observe that we can 
perform the change of variables r — (y/1 — M/ir)f, ip = — M/n)(j> and observe that close to r = the geometry 
of a slice at constant time is given by 

df 2 + r 2 d(p, 

with < <p < 2ir(y/l — M/tt). This corresponds to a conical singularity of deficit y/1 — M/tt. From (p~2|) one obtains 
that < M/tt < 1, which is the physically admissible interval^ (M/tt > 1 corresponds to sources of negative 
energy-density). In fact, from the definition of 7 one gets 

1 + T 2 - = (2M*+-7) 2 , 

and therefore 

(l + 7 2 )>(2^ + - 7 ) 2 . 

It is also clear form the definition that 2[il + > 7, and therefore we obtain what we need, that is, 



\f\ + 7 2 + 7 = x > 2[il + . 

The value of B in Eq.([T]) can be readily computed using the Hamiltonian formulation. This is done in the appendix 
]where we find that, on the constraint surface, the action is given by, 



I = -H I \fhd 2 o- ( ^E 2 d 3 x+ [ *F AA + nA+l- [ (N*f lP -l). (13) 



Now, 



where 



*FAA = E+I A = E+I F = El/ ^ = £ 2 2tt/3_ / rdr = n[3-, (14) 

dM JdM JM JM JO 



+ 



/3_ = 2-nl-y/X -M/tt (15) 

is the inverse temperature associated to the cosmological horizon of the (— )-region. We obtain a final expression for 
B given by 

B = -I + I Back =TT(l3+-f3_) (16) 

where [3+ = 2ttI+ is the temperature of horizon of the (+)-region. 

Recalling that e~ B stands for the tree level semiclassical approximation to the probability of nucleating a string 
through thermal fluctuations we have found out the final expression to be evaluated. 

Assuming that q > we notice, from Eq. (|5|6[l . that /_ > 1+ thus the cosmological constant is reduced in the region 
r > R when the string materializes. After the nucleation the system may be analyzed classically. In the Lorentzian 

sector, the thermalon is unstable, and may collapse both towards r = r or towards r = 0. In the first scenario, 

the (+)-rcgion grows to cover the entire space so that the initial situation is recovered. In the second one, the string 
collapses into a point particle [3] and the (-)-region fills the entire spacetime. This is analogous to what happens in 
four-dimensional spacetime where the final geometry is a Schwarzschild-de Sitter geometry with a smaller cosmological 
constant. More details on that process in four dimensions may be found in [3j|. 
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III. THERMALON IN FIXED CURVED BACKGROUND. ANALYSIS OF PREFACTOR. 



We will start by analyzing the operator O in {2} for the case in where the string back reaction in the gravitational 
field is neglected, i.e., in the limit G — ► 0. The first consequence of this is that separation into ± sectors becomes 
irrelevant since neither the geometry nor the cosmological constant jump across the string, thus l + = I- = I. 

Unlike the geometry the electromagnetic field jump across the -charged- string in concordance with the Gauss 
constraint ([5]) . However in this case the electromagnetic field does not incorporate degrees of freedom since it is fixed 
between r = and the position of the string and determined by ([5|) behind the string. Therefore the only degrees of 
freedom are those of the string. 

Subtracting the background to the action ([3]) we find 



I — Iback — Ithe 



V9 d3 



i\f~hd 2 o, 



(17) 



where e = E\ — E 2 _. The first integral is between the radius of the string R(t, <fi) and the cosmological horizon /. In 
the second term hij is the induced metric on the string. For simplicity we choose the coordinates a = (t, </>) on the 
string. The induced metric is 



dh 2 



f 2 (R) 



R 2 



P(R) 



dt 2 



R 2 



R' 



f 2 (R) 



„ RR' 

2——dtd<b. 

f 2 (R) 



Integrating the first term in (|17[) in the radial variable we find, 



dtd<fi 



R 2 



Vh 



(18) 



(19) 



We now insert the de Sitter metric ([7]) with cosmological constant 1/Z 2 and the induced metric (JTSJ) in (fT9| . The 
resulting action is function of the position of the string only, R(r,(f>). In order to compute the number of negative 
eigenvalues of the operator O we will expand the around the classical -static- solution of radius Rq as 

R = Ro + A(t, <j>) (20) 

up to second order in A such that 

I = I + Ii+I 2 +O{A 3 ). (21) 
The first term, Iq = —B, was computed in the previous section and it is recovered here. The second term, 

h = -A (|/(i? )iW 2 + Kl 2 - 2iZ§)) , (22) 

vanishes when the equations of motion are used as expected. It is easy to see that I\ is equivalent to Eq. (|ll[) in the 
limit G -> 0. 



The operator O can be read from, 



I 2 = -A 2 (l + ^(2^^ 



which using Eq. 



2 ' P 
may be rewritten as 

h = 



. Rn 



2/i 

7^ 



- A 



(A* - R 2 A 2 - f 2 A rz ). 



[iRc 



-A n 



Rof 



2 a /2n 



(23) 



(24) 



Making the change of variable t — £-t, with < t < 2ir and j3 = 2ttI the inverse of temperature of the cosmological 
horizon, the operator ([2]) is given by 



O 



4^ 

PRo 



dl + fd 2 



(25) 
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The eigenvalues of O are proportional to the positive constant 4/i// 3 i?o. As we are only interested in the sign of the 
eigenvalues, we will only consider the ones associated with the operator in square brackets of (f2"5| . They are 

X nm - -1 + n 2 ^ + m 2 f = -l + m 2 + ^(n 2 - m 2 ) . (26) 

where n, m G Z. 

Recalling that Rq < I 2 , it is direct from Eq. (j2"6")) to notice that there are five negative eigenvalues: Aoo, A±io and 
Ao±i- Furthermore, there are four null eigenvalues given by A±i±i = 0. The presence of these null eigenvectors is due 
to the 50(4) symmetry of the background geometry. Indeed, from the six symmetry generators, only four rotate the 
string. The remanning two rotate it on itself, corresponding to reparamctcrizations of the worldshcct. In appendix lAl 
we show that it is also possible to count the number of negative eigenvalues making use these symmetries. It is worth 
noticing that this result is in complete agreement with the mathematical analysis carried out in [(| where the number 
of negative eigenvalues of the operator associated with area of a two-torus immersed on S 3 is, at least, five. 

If we accept an arbitrary value of (5 then we may reduce the number of negative eigenvalues to only three. In fact, 
it is easy to see that if 

[3 < 2ttR (27) 

the only eigenvalues which are always negative are Aoo and Ao±i- This means that if the temperature is high enough, 
i.e., if one puts the system in contact with a heat reservoir with a temperature satisfying (|27p . then one may reduce, 
in the best case, the number of negative eigenvalues to three. This, however, amounts to admit a conical singularity 
at the cosmological horizon. 



IV. THERMALON INTERACTING WITH GRAVITY. ANALYSIS OF THE PREFACTOR 

Now we let gravity interact dynamically with the string. We are interested in computing the one-loop correction 
to the semiclassical approximation of the path integral corresponding to the following Hamiltonian action, 

1=1 •://,,.-'' • .i,,/"' A H \'H, A,C"i./ :! r • f d 2 xp^ + nA. (28) 

Here gij, tt^ are the canonical variables describing the gravitational field, Aij, P u are the ones corresponding to the 
electromagnetic fields, while the pair y^, p^ correspond to the string degrees of freedom. TV 1 , Tt ± are the gravitational 
constraints, and G l the electromagnetic Gauss constraints. N % , N 1 - ,Xi are Lagrange multipliers. The last term is 
A/AGh, A being the area of the cosmological horizon, and it is required for ensuring invariance under change of 
coordinates (We use the same conventions adopted in pj|). 

The action (j2"8|) contains gauge degrees of freedom. One way to perform the path integral is to fix the gauge and 
work in the reduced phase space. This amounts to include a delta function on a set of gauge fixing conditions, 

x M = o 

and to complement the integration measure with the determinant of the operator x^Ksee, for example, jll|). 
Here we run into difficulties. In the gauge we have chosen (see appendix [Bj, the determinant in question turns out to 
be untractable. On the other hand, this gauge has the advantage of allowing us to completely solve the constraints 
and write the reduced action explicitly. We, therefore, will work with the reduced action and will assume that the 
new measure will not affect the number of negative eigenvalues we are counting. 
The metric may be written in an ADM form as follows, 

ds 2 = N 2 dt 2 + dp 2 + 7 2 (d0 + N^dt) 2 . (29) 

Here we have chosen the following gauge fixing conditions 

9pp = 1 a P 4> = 0. (30) 

To completely fix the diffcomorphism invariance we still have to choose a third condition. We require 7r^* = 0. These 
gauge fixing conditions were first used to completely solve the constraints in [l2j]. We are going to gauge fix the 
electromagnetic field as well by requiring A^ p = 0. 
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The constraints TV = Ji 1 - = 0, whose explicit form is included in appendix [Bj along with the gauge conditions form 
a set of second class constraints that may be completely solved outside the strings, in terms of four functions of (t, (f>). 
At each side of the string, we get 



r 2 cos 



P~ Po 
I 



r+ sin 



P- Po 
I 



Ik 

1 d (r- 
- — — arctan — tan 
Kocp \r + 



(31) 



P~ Po 

I 



n. 



The four functions r+(t, 0), r_(t, <fi), po{t, <ft), n(i, (f) are arbitrary. However, in general, they will be different at each 
side of the string. The way they jump across it is determined by the constraints, as we shall see below. 

The conservation of the gauge conditions (f30|) in time, — g^p = g pp = 0, give rise to a set of consistency 
conditions. 







N p 



■ 1 2 N*+2k-^N 



1 



(32) 



= 



These equations restrict the possible values of the Lagrange multipliers N, N l . The most general solution may be 
given in terms of four arbitrary functions of (t,<ft): u(<p,t), v(4>,t), w(4>,t) and n p ((j),t), 



N p = 



1 






MCOS 
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w - 


1 






n" 





Po 



sin 



Po 



P~ Po 

I 



I 

2r 2 _r_w sin 



v r i sin 



P~ Po 



P~ Po 

I 



— r cos 



P - Po 



P- Po 

I 



I 



(33) 

(34) 
(35) 



where again the fields may jump across the string. 

To compute the operator associated with the second variations of the action, we consider the most general variation 
of the metric and string position respecting the constraints above. 

The variational principle we use is adapted for fixing 7, N l and N at the boundary p = 0. We fix 7 = 0, = 0, 
N p = and N = 1, therefore the metric in the (+)-region has 



T(+) 



= r.(+W ^ 



N 9 =0 
iv (+) ' 



N, 



(+) 



(+)«, (+), 



+ 



(36) 
(37) 
(38) 



The fields r_W and 

remain arbitrary and we have to sum over them in the path integral. 
The metric functions in the (— )-region have six functions determined by the three continuity conditions for 7, N 
and and the three jumping condition (|B14|) . (|B18j) . (|B19j) for its radial derivatives across the string. This is, in 
general, a very difficult task. However, we only need to consider the second order variations of the action, therefore 
we make an expansion of all functions around the thermalon solution. For the position of the string we write, 



P 

P 
-1 
P 



+ Xp 



(i) 



(39) 
(40) 
(41) 



In terms of this parameter we write the variation of the remaining fields to order A 2 , 
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rJ +) - Z_l 4- XrJ 1+) + X 2 rJ 2+) +0(\ 3 ) 


(42) 


- Av (1+) + AV 2+) + OCA 3 ) 


(43) 


rJ-t = XrJ 1 ^ + X 2 rJ 2 ^ +C(A 3 ), 


(44) 


r+ <-) = >/MZ_ + XrJ 1 ^ + A 2 rJ 2 ^ + 0(A 3 ), 


(45) 


= CVM^ + Xu^ +X 2 u^ 2 - ) +0{X 3 ), 


(46) 


«(-) = A^ + AV^+OtA 3 ), 


(47) 


u,(-) = Aw (1 ~ } + aV 2 ~> +0(A 3 ), 


(48) 


Po 1 "' = p + Ap (1 - ) +A 2 p ( 2 -)+0(A 3 ). 


(49) 


der by order in A. We end up with five arbitrary functions, p^ 


p(D } p'd), r+ d+) 



and in terms of which all the other fields are determined. Each of these five functions must be integrated in the 
path integral. 

Finally, the action evaluated on the Maxwell equations (IB6[) is given by 



-Ei 



N(+)~t( + )dpdtd4 - Ei 

CP+ + 



N ( -)1(-)dpdtd(/> - n / \jN 2 p' 2 - N 2 j 2 - p>y 2 dtd<f> 



+ i/ K-) dtd< t> + E + J o ++ N {+)1{+) dpdtd<j } +^ j '(JV (+) 7(+), P - l)dtdcj>, 



(50) 



where the third term is evaluated on p = p, the fourth term is evaluated on p = p (the cosmological horizon), the 

last term is evaluated at p = and p++ = ttI + /2 is where 7+ is maximum. The fourth term, on shell, is equal to 
A/ 4 where A is the area of the cosmological horizon. Replacing the functions as series of A in ([50]) we can express 
the action in the following form, 



I = I + Xh + A 2 / 2 , 



where the term I2 can written as, 

h 



~ (p (1 \a + bd 2 + cdl)p {1) +UxU + VyV^j 



(51) 



(52) 



Here the parameters a,b,c,x and y are functions of l + . Z_ and p. The functions U and V are linear combination of 
pW, r+( 1+ ) and v^+l 

The operator O can be written in this case as 



O 



f a + bd 2 + cdl 




The eigenvalues of this operator are 



Xnr 



Ai = x, A 2 = y 
2vr N 
J 



(53) 



(54) 



b ( — - I n 2 — cm 2 



where n, m G Z. 

We study numerically the quantities a, b, c, x and y for different values of the parameters. The result of those 
numerical computations yields the minimum number of negative eigenvalues is twelve. This result proves that the 
presence of gravity increases the number of negative eigenvalues. 



V. CONCLUSIONS 



In this work we have analyzed the thermalon process in three spacetime dimensions, centered in the problem of the 
number of negative eigenvalues of the operator O that arises from second order perturbations around the classical 
solution. 
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When gravitational back reaction is neglected, the operator has five negative eigenvalues. We point out that if the 
fixed geometry represents a de Sitter geometry in thermal contact with a reservoir with high enough temperature, 
then the number of eigenvalues may be reduced to three. This, however, amounts to accept a conical singularity in 
the background geometry. 

Since the presence of more than one negative eigenvalue seems to rule out the interpretation of the thermalon as 
the indication of metastability also in this work we explored the influence of gravity expecting that this would reduce 
the number of negative eigenvalues. The result was the opposite and the presence of gravity seems to increase the 
number of negative eigenvalues. 

In this work we neglected the influence of the measure in the path integral. Although the correct treatment of 
the measure may reduce the number of negative eigenvalues, as it does for the conformal mode of gravity in [l3j . in 
our case it is unlikely that this could account to eliminate a plethora of negative eigenvalues. Therefore, we argue 
that the inclusion of the gravitation backreaction not only does not reduce the number of negative eigenvalues but 
most probably it actually increases them, spoiling the standard interpretation, a la Coleman, of the thermalon in the 
context of thermal decay rates. 
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APPENDIX A: ANALYSIS FROM THE SYMMETRIES 



In this appendix the result in Eq. (|26p will be reobtained starting form the symmetries of the problem. This is 
in analogy of what was done by Garriga in 14|. We will study the fluctuations of the string around the thermalon 
solution r = Rq. The operator that represents the second variations in the action must be covariant on the string 
worldsheet. This implies that it must be proportional to —V 2 + C, where —V 2 is the laplacian on the worldsheet and 
C is a constant. Now, 



1 



1 



V - —d 2 + — & 



Rq 



(Al) 



We define a dimensionless coordinate r 



2^, so that t S 



0, 2tt], and the operator is proportional to 



p(l 2 ' 



c. 



(A2) 



The eigenvalues in the space of 27r-periodic functions in <f> and r are, therefore, proportional to 



(2^) 2 n 2 
~1W 



R l 



C, 



(A3) 



where n, m £ Z. In order to determine C we note first that this constant does not depend on f3. This is because 
the global parameter /? appears only through the redefinition of the time coordinate, and therefore the Lagrangian 
density, which is a local quantity, cannot depend on it. Hence, C may be computed for some particular value of fj 
where the answer is known. In the case in which [3 — 2nl, the geometry is a 3-sphere (globally Euclidean de Sitter) and 
one knows that 4 zero modes must show up. This is because the 3-sphere has six symmetries, two of which leave the 
position of the thermalon, r = Rq, invariant (representing reparameterizations of the worldsheet of the string). The 
other four should give rise to these zero eigenvalues. Note that the corresponding killing vectors must have periodicity 
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2ir, for they correspond to rotations, and therefore they must have n = m = 1. This are precisely the four modes 
we need, and therefore 



which implies 



C ""W (A5) 



As a double check let us take another limit where we know the number of 0-modes: flat spacetime at finite temperature. 
This is 3 dimensional flat space with the Euclidean time coordinate identified with periodicity (3 which correspond to 
the inverse temperature. The complete six-dimensional group of rotations and translation is broken to translations plus 
the rotation in the plane orthogonal to t. The other two rotations do not respect the identification (do not commute 
with dt). From the four remaining symmetries, two are, again, just reparameterizations of the string worldshcet, and 
therefore we are only left with two spacial translations, which correspond to killing vectors independent of t that have 
periodicity 2tt in space. We expect that in this case the eigenvalues with n = 0,m = ±f vanish. In fact, 



lim A 0± i = \ + lim C = 0. (A6) 
R 



Now plugging (|A5|) back in (|A2|) we obtain precisely our result (f2"5)) of Section . 



APPENDIX B: THE HAMILTONIAN FORMULATION OF THE THERMALON 

The Hamiltonian action is given by (|28p . On the constraint surface, the canonical action is simply 

1 = J p^ + ttA. (BI) 
The canonical momentum p^ conjugate to the string is is given by 

p, = -vVhg^(h 00 y u + h°Wn - qA^y'r (B2) 
Replacing Eq. (|B2|) in Eq. flBljl this action may be written as, 

I=-H [ d 2 aVh -q [ A + it A. (B3) 

The action principle (|28[) must be supplemented with information on the fields which are going to be fixed on the 
boundary. In the present case the boundary is defined to be at r = 0, mimicking what was done in ||. When varying 
the action one gets the following boundary terms at p = 

-~ / (N6 %P ~N P 6 1 )+ [ N^e. 

z JdM JdM 

We see that the action is ready to fix 7 and 7 P and 7rf on the boundary. We are interested in fixing N = 1 instead 
of 7 P , so that the inverse temperature = (t% — t\) will be fixed, as it should be in the canonical ensemble. To 
accomplish that we may add to the action the boundary term 

~/ (N% P -K), 

L JdM 

where K is an arbitrary constant. In order to compute K we consider the situation when there is no string and we 
have a particle singularity at p = r = 0. In that case the boundary term on shell is given (3(1 — M — K), which shows 
that K = 1. In fact, the complete action on-shcll is 



(3(—M) - it A. 
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Recalling that the action in the canonical ensemble represents the free energy, and ttA the entropy, we find the internal 
energy is —M. This is in agreement with what was done in [l5j . The mass will be +M when the boundary is set at 
the cosmological horizon. 

After taking the boundary terms into account, the action takes the form 



1= -/i f cfaVh-q [ A + tiA+\ [ (AT 7 , p - 1). 

JS JT. * JdM 



The electromagnetic constraint in the action requires P lJ to be fixed on the boundary. In our case, we want 
P*i = E + £ l; ! at the boundary, where e 4 - 7 is the Levi-Civita tensor on the constant t slice. The constraints require that 
E+ is independent of <p. Furthermore, we want to describe a time-independent cosmological constant and therefore 
we choose E+ to be a constant. Gauss' constraint G % will completely determine the value of the field E everywhere. 
We hence find that on the constraint surface we are allowed to use the Maxwell equations d*F = — *J. We may then 
write, on the constraint surface, 



q I A= *JAA = 

'£ J M 



d*F A A = *FAF- *F A A, 

M J M JdM 



where 8A4 is the boundary of M. at r — 0. The action is 

/ = —(J, I VhcPcr — f *F A F + f *F A A + nA + i f (N^ p -1), 
Jt, J m JdM 2 J dM 

which may also be written as, 

I = -H I \fhd 2 a- I ^g~E 2 d 3 x + I *F AA + tiA+ \ [ (iV 7 , p - 1). 

JT, J M JdM 2 Jg M 



(B4) 



(B5) 



(B6) 



We finish this section by computing the form of the gravitational degrees of freedom left after solving the constraints. 
The Hamiltonian constraints contain a pure gravitational piece and a matter term, H = Ji Grav +7i Ma *. In the present 
case 



n 
n 



Grav 
P 

Grav 



-2(< p + { 7 - 2 ^}, ) , 
-2ttP „. 



The matter part is given by 



H 



Mat 



y/gT ±± and U\ 



Mat 



_Li, 



(B7) 

(B8) 
(B9) 

(BIO) 



where ^fg = 7, Tj_j_ = T^n^n^, T±i = T^iti^ and = (N, 0,0,0) is the unitary timelike vector normal to the 
t =constant surfaces. The string may be parameterized by p = R(t, (f>). The induced geometry on it is 



da 2 = (N 2 + R 2 + ^N^dt 2 + 2 (Aft' + 7 2 N^dtdcj) + (R' 2 + 7 2 )# 2 . 
The energy momentum tensor produced by the string and the electromagnetic field is given by 

1 



T" v = ^E 2 g^ + n 



(3) da a da h 



(Bll) 



(B12) 



'9 



Now, using (|B10j) - (|B12j) . we may compute the matter part of the constraints, 

Vh 



H 



Grav 



7 3 iV 



E 2 
2N 2 



J ^S^(x-y)d 2 a 

H A J at = 1 N[i [ y^(h tt R + h^R')S^(x-y)d 2 a 
' J Nj 



(B13) 
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To find the jumping conditions we integrate the constraints and consistency conditions across the string. We obtain 

2 A R' 2 + 7 2 , , 

~k p = M ~ 

2 

2A7r" p = ^L(R-R'N't') (B15) 



. (RR' + A^ 7 2 ) 



2A7r" = - 7 V ^= ; - (B16) 



For the TV's, we note that the momentum 7r p in terms of the metric is given by 



*** = ^V ( B17 ) 



3 



27V 



Replacing (|B17|) in (|B16|) we obtain 



AN*, P =-^ {kR ' + ^ 2) . (B18) 
7 V/i 



In the same way, (|32[) may be integrated around the position of the string. This gives the discontinuity of the 
derivative of N, 

-AN, p =-N 2 1( i-^(y/hh tt ) (B19) 

Note that in the expressions (|B14[) . (|B18|) and (|B19[) the right hand side is evaluated at the position of the string. 
This is not a problem, for N, 7 and are continuous functions. 
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